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Abstract 

We derive an exact representation of the massless Schwinger model on the lattice in terms of 
dual variables which are configurations of loops, dimers and plaquette occupation numbers. When 
expressed with the dual variables the partition sum has only real and positive terms also when a 
chemical potential or a topological term are added - situations where the conventional representation 
has a complex action problem. The dual representation allows for Monte Carlo simulations without 
restrictions on the values of the chemical potential or the vacuum angle. 



1 Introduction 


Since its initial formulation three decades ago lattice QCD has developed into a reliable quantitative 
tool for studying many low energy phenomena in QCD. However, one important issue where the lattice 
approach has essentially failed so far is its application to QCD at finite density. The reason is that 
at finite chemical potential // the fermion determinant becomes complex and cannot be used as a 
probability weight in a Monte Carlo simulation. This is known as the "complex action problem” or "sign 
problem”. The complex action problem is not specific for QCD or for theories with fermions, but is a 
generic feature of many quantum field theories at finite density both in the lattice and the continuum 
formulation. Furthermore, not only finite chemical potential gives rise to a complex action problem, but 
also the addition of a topological term. 

For overcoming the complex action problem in lattice simulations several approaches have been 
discussed, such as reweighting, various series expansions, the reformulation with alternative variables 
and simulations based on stochastic differential equations. Reviews at the yearly lattice conferences [T| 
give an overview about the progress in recent years. 

Probably the most elegant approach to the complex action problem is to rewrite a lattice field 
theory in terms of new degrees of freedom, often referred to as "dual variables”, such that in the new 
representation the partition sum has only real and positive contributions. The dual variables usually turn 
out to be loops for matter fields and surfaces for gauge fields. The key problem of this approach is that 
there is no general recipe for finding a dual representation of a given theory. In particular non-abelian 
gauge fields are difficult and so far no convincing approach has emerged. The situation is different for 
theories with abelian gauge fields where for the case of their interaction with bosonic matter several 
interesting results based on dualization [2l [3] were presented in recent years. 

Coupling fermions to the gauge fields brings in additional difficulties with extra signs for the matter 
loops coming from the Grassmann nature of the fermions and from traces over gamma-matrices. In 
2 dimensions, however, the understanding of the traces over gamma-matrices as they appear in dual 
fermion loops is much more advanced |l4|, such that the Schwinger model, i.e., QED in 2 dimensions, 
seems to be an interesting candidate for further developing the dual approach for theories of gauge fields 
interacting with fermions. Indeed, some results for at least partial dualizations of the lattice Schwinger 
model were discussed in the literature [5l[6l|7], but no convincing complete solution of the sign problem 
was presented so far. 

In this article we show that for the lattice Schwinger model with massless staggered fermions the 
complex action problem can be solved completely by exactly rewriting the model to a dual representation. 
We consider both types of complex action problems, coming from finite density and the addition of a 
topological term. The dual variables are closed oriented loops for the fermions, with the gauge fields being 
represented by integer valued plaquette occupation numbers. The constraints for these dual degrees of 
freedom again allow for the interpretation of the dual gauge field variables as surfaces in 2 dimensions. 
We show that in the dual form the partition sum has only real and positive contributions. Furthermore 
we discuss the dual formulation of various observables and briefly address possible update strategies. 

We stress that the dual representation discussed here not only provides an interesting step towards 
solving complex action problems for theories of gauge fields interacting with fermions, but that it will 
also be useful for addressing some of the open questions concerning the phase structure and universality 
class of the massless lattice Schwinger model with staggered fermions |E[9|. 
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2 Conventional form of the one flavor model with topological term 


We begin with the case of the massless one-flavor model with a topological term. In the conventional 
representation the corresponding partition sum is given by 


The dynamical degrees of freedom are U(l)-valued link variables f7t/(n) = exp(z^,y(n)), Ay{n) G [—tt, vr]. 
Here n = (ni,n 2 ) denotes the sites of a 2-dimensional Ns x Nt lattice and u = 1^2 runs over the 
Euclidean space [v = 1) and time {v = 2) directions. The second set of degrees of freedom are the 
one-component (we work with staggered fermions) Grassmann valued fermion variables 'ip{n) and ip{n). 
All boundary conditions are periodic, except for the temporal boundary conditions of the fermions which 
have to be chosen anti-periodic. By V = NsNt we denote the total number of lattice points. 

The measures in the path integral are the product measures of U(l) Haar-measures for the gauge 
fields and Grassmann measures for the fermions: 


n 


dAy{n) 
- 2'n 


n/ dip{n) d'ip{n) 


( 2 ) 


For the gauge action we use the Wilson form (the constant term was dropped for simplicity), 


Sg[U] = -pJ2^eUp{n) = -^^[[/^(n) + Up{n)-^] 


(3) 


The plaquette variables Up{n) = Ui{n) U 2 {n + 1) Ui{n + 2)“^ U 2 {n)~^ are the usual products of the 
link variables. We remark, that in two dimensions all plaquettes Up{n) can be labelled by the coordinate 
n of their lower left corner. For the topological charge Q[U] we use the field theoretical definition 


= J^T.lUp{n) - Up{n)-^]. (4) 

We stress that the field theoretical definition Q approaches the continuum form Q = \/2'k f d‘^x Fi 2 {x) 
only in the continuum limit. At finite lattice spacing Q[U] is not an exact integer and 6 is not an angle 
and consequently requires renormalization. In a recent publication f3l| it was shown non-perturbatively 
in a lattice simulation of scalar QED 2 how a suitable continuum limit can be set up, such that Q[U] 
as defined in converges to the continuum form and the 0-dependence of observables becomes 27r- 
periodic, i.e., 0 indeed becomes an angle. 

The action for massless staggered fermions is given by 




+ “ lv{n)Uy{n) ^ ii{n + i))'4^{n) , 

n,z/ 


(5) 


where the staggered sign function ^ij{n) is defined as 


7 i(n) = 1 , 72 (n) = (-1)"^ . 


( 6 ) 


Putting things together we arrive at the following form of the partition function 






( 7 ) 
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where we have introduced the abbreviations 


V 


^ _ A 

2 ~ 4 ^ 


_ ^ e 

11 = —I- 

' 2 47r 


( 8 ) 


It is obvious from Q that the conventional representation has a complex action problem for non-zero 
values of the 9 angle. Then rj ^ rj and the Boltzmann factor acquires a phase and is therefore not 
suitable as a probability weight for importance sampling. We remark, that for the physically interesting 
region where we want to construct the continuum limit (/3 —)• oo) both parameters r] and fj are positive. 
We will assume this positivity from now on, i.e., we restrict ourselves to /3 > 1/2 with 9 G [—vr, vr]. 

For later use we redefine the link variables Uu{n) by absorbing the staggered sign function 7 j/(n): 


Ul{n) = -iu{n) Uu{n) . (9) 

This transformation obviously removes the staggered signs from the fermion action (|^. The plaquette 
simply changes sign, 

Upin) = U[{n)U' 2 {n + i)U'p{n + 2)-^U!2{n)-^ = -Ui{n)U 2 {n + l)Ui{n + 2)-^U2{n)-^ = - Up{n). 

( 10 ) 

Due to the fact that the values ±1 of the staggered signs are elements of the gauge group U(l) the 
path integral measure for the gauge fields remains invariant {^T>\U'] = /D[[/]) and we obtain for the 
partition sum (for notational convenience we drop the primes on the transformed link variables) 

z = jv[U] e-^^n^vin) Z^[U] = jv[U] Z^[U] , (11) 


where we have defined the fermionic partition sum Z^[U] in a background gauge field U as 


Z^U] 




Uu(n) ^ 


( 12 ) 


With the transformation ([^ we have removed the staggered signs and traded them for an extra minus 
sign for the gauge action (and the topological charge term). 

We remark that the transformation can be applied in arbitrary numbers of dimensions and for all 
theories where —1 is in the gauge group (U(l), SU(2N) et cetera). It always has the same effect of 
removing the staggered signs in the fermion action and flipping the sign of the Wilson gauge action. 


3 Integrating over the fermion fields 


The next step is to integrate out the fermions in the fermionic partition sum Z^[U] defined in (12). 
For that purpose the Boltzmann factor in Z^[U] is expanded, such that the Grassmann integral can 
be saturated with the terms from the expansion. The integral is saturated when at each site n both 
Grassmann variables ^/(n) and ili{n) appear in the integrand as factors generated from the expansion. 
In other words, the only non-vanishing Grassmann integral is 


il){n)'4>{n) = / dV’(n)dV'(n) ^/(n)V'(n) = 1 


(13) 
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We proceed by rewriting the sum in the exponent of the Boltzmann factor of (12) into a product and 
expand the individual exponentials 


Zf[U] = 


-5 C/i/(n)'0("•) ■0("+£') p\Ui,{n) ^' 4 ){n+i/)' 4 ){n) 


V'] e 2 '"''VV rvv vv-Ig 2 


n,u 


V'] n E (- 

n,u k„{n)=0 


U fr?) - ^ — 


_i: I 

k^{n)=0 


Ujn)-^ - 


(14) 

kv{n) 


V^(n+i>) ip{n) 


Since the Grassmann variables are nilpotent, the power series for the exponentials terminate after the 
linear term. Thus the summation indices ki,{n) and ky{n) which we use for the forward and backward 
hopping terms run only over the values 0 and 1, and we will refer to them as "activation indices” for the 
corresponding terms. The terms that are activated, i.e., those with ky{n), k^{n) = 1, bring as factors 
the corresponding link variables Uy{n) for forward hops and Uy{n)~^ for backward hops, as well as trivial 
factors of 1/2 (which, as we will see, combine to an overall factor of (1/2)'^). In addition all activated 
forward hops contribute a factor of —1. 

The simplest way for saturating the Grassmann integrals on two neighboring sites n and n + i> 
is to activate the forward and the backward hopping terms that connect the two sites, i.e., one sets 
ky{n) = ky{n) = 1. Such a term is referred to as "dimer” and the corresponding contribution is 


Uy{n) Uu{n) ^ 
2 


Jdip{n+C')d'ijj{n+h)dijj{n)dip{n) r/(n)r/(n+z>)'i/(n+i>)?/(n) 


1 

4 ■ 


(15) 


The link variables cancel and the overall minus sign which comes from the forward hop is compensated by 
a minus sign from the interchanges of Grassmann variables when the last factor ?/(n) of the integrand is 
commutated to the beginning of the integrand in order to have the Grassmann variables in the canonical 
ordering of the integral ( fls] ). Thus, up to a factor of 1/4, which we will absorb in an overall constant 
for the partition sum Z^[U], a dimer saturates the Grassmann integrals trivially. 

The second type of contributions that saturate the Grassmann integral are oriented closed loops, 
where the activation indices k^{n) and ku{n) along the contour of a loop are equal to 1. Obviously 
a hopping term in the loop provides a 'ijj{n) and the next hopping term in the loop the corresponding 
^p{n) to saturate the Grassmann integrals for all sites in the loop. We obtain an overall minus sign from 
commuting the last 'ip{n) to the very beginning of the integrand. Thus each loop comes with an overall 
minus sign. In addition every forward hop in the loop comes with a minus sign and, since for a closed 
loop the number of forward hops is half of the total number of hops in the loop, we can write this sign as 
(_l)i(0/2, where L{1) denotes the length of the loop 1. This counting of the signs from forward hops is 
also correct for loops that close by winding around the boundary if we choose the number of sites in each 
direction to be a multiple of 4 (otherwise trivial corrections have to be taken into account). For loops 
that close around the compact time direction we have an additional minus sign from the anti-periodic 
temporal boundary conditions of the fermions, but we will show below that the total temporal winding 
number of all loops must be even and this sign will turn out to be irrelevant. Furthermore, each hop 
comes with a factor of 1/2 and the corresponding link variables Uii{n) for forward hops and Uy{n)~^ for 
backward hops. Thus the overall contribution of a loop I is: 

(_l)m0(_i)|i(0 (16) 
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Figure 1: Example of an admissible configuration on an 8 x 8 lattice. In the Ihs. plot we show only 
the fermion contribution, i.e., the filling of the lattice with dimers (the double lines on the links) 
and with oriented loops. In the rhs. plot we also show the plaquette occupation numbers needed to 
compensate the link flux introduced by the loops. 


As before L{1) is the length of the loop I, i.e., the number of links in the loop and W{1) £ is the 
winding number around the compact time direction. The product is over all links {n,iy) in the loop I, 
and Su{n) takes into account the orientation the link is run through in the loop, with s^{n) = I for 
forward hopping and s^{n) = —1 for backward hopping. 

For completely saturating all Grassmann integrals we can now combine loops and dimers such that 
each site is either run through by a loop or is the endpoint of a dimer. Note that due to the absence of 
mass terms, there are no monomer terms and a non-vanishing contribution to Z^[U] must saturate the 
Grassmann integrals at each site with either a dimer or a loop. Thus Z^[U] is a sum over the set {l,d} 
of all configurations of loops and dimers: 



where denotes the total number of loops in a configuration. The factors 1/2 from the hops in the 
dimers and in the loops combine into a trivial overall factor of (1/2)^. 

In the Ihs. plot of Fig. [T] we show a legitimate configuration of loops and dimers that contributes to 
the fermionic partition sum dlT] ): Each site is either run through by an oriented loop or is the endpoint 
of a dimer. Note that this constraint restricts the possible shapes of the loops, and, e.g., a 2 x 2 loop 
is not possible since the single site inside the loop cannot be saturated by a dimer. In more than two 
dimensions or in the presence of monomers from mass terms a 2 x 2 loop would be possible. 


5 












































































































4 Integrating over the gauge fields 


After having found the representation of the fermionic partition sum Z^[U] in terms of dimer and loop 
configurations, we now can integrate out the gauge links. The gauge links come from two different 
sources: Each loop in is dressed with the link variables Uy{n) sitting on its contour, and from 

,-^Upin)-^ 


in (11) we obtain further factors of 


the expansion of the Boltzmann factors and 

link variables. Collecting all contributing factors each link variable Uu{n) appears with some power 
jv{n) G Z. Using the well known representation of the Kronecker-Delta, l/27r = 5jfl for 

j G Z, we find the master formula for the gauge links. 


TT = TT' 


(18) 


The formula implies that for a non-vanishing contribution to the path integral only terms where all links 
Uu{n) appear with exponent ju{n) = 0 are admissible, i.e., all link fluxes cancel. Thus every link variable 
Uu{n) that is generated by a loop in Z^[U] has to be compensated by a factor Uy{n)~^ from expanding 
the gauge action terms. In other words, we need to occupy the plaquettes inside the loops by powers o’ 
Up{n) generated from the expansion of ^-vUp{n) i jg jiijstrated in the rhs. plot of Fig.jl 

where we indicate the occupied plaquettes by their occupation numbers, as well as by arrows to indicate 
how the link fluxes cancel. The arrows of a plaquette at n are oriented in the mathematically positive 
(negative) sense for p(n) > 0 {p{n) < 0). No arrow indicates p(n) = 0. The figure also illustrates that 
for links inside a loop the link flux introduced by adjacent plaquettes cancels. 

At this point we now also see that the sign (—from loops with non-vanishing temporal 
winding is irrelevant for admissible configurations: In order to saturate all gauge links we need to have an 
equal number of forward and backward winding loops, since adding plaquettes only shifts the unsaturated 
links of a single winding loop. Thus the total number of forward and backward winding loops is always 
even and from now on we can drop the factor (— 

For a convenient algebraic treatment of the gauge integration we use the following representation of 
the Boltzmann factor for the plaquette Up{n), 


^-'nUp{n) ^-r]Up{n} ^ 


p(n)ez ^ 


(19) 


where I\p(n)\{‘^V^) denotes the modified Bessel functions, and p{n) G Z is referred to as the plaquette 


occupation number for the plaquette at site n. The rather elementary proof of (19) is given in the 
appendix. Note that for our restriction to positive rj and p the argument of the Bessel functions is real 
and positive, implying that also the values (2-y/r)f)) are real and positive. 


Putting things together, we obtain the following result for the full partition sum. 


'.2 


E(-i) 


Nl+Np + P 


{14,p} 


(Vf 


p{n) 


where we have introduced the total plaquette occupation number 

n 


( 20 ) 


( 21 ) 


The sum J2{i,d,p} '>1 (20) runs over all admissible configurations of loops, dimers and plaquette occupation 
numbers. The admis^le configurations have to obey (see the rhs. of Fig. [^for an example): 
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Figure 2: Example for the factorization of a configuration with nested loops. The configuration is 
decomposed into two factors, where in each factor one of the two loops is replaced by a closed chain 
of dimers. 


• All sites of the lattice have to be either run through by a loop or be the endpoint of a dimer. 

• At all links the fluxes introduced by the fermion loops and from occupied plaquettes must cancel. 


Some remarks are in order here: For non-vanishing vacuum angle 6 ^ Owe have rj ^ rj, and the factor 


p{n) 


in (20) gives a different weight to positive and negative plaquette occupation numbers. 
However, except for the explicit sign (—+ + 2 Si weight factor is real and positive, i.e., 

Yin ^ ^ ®'+' Thus, if one can show that the sign (— 1 )^^+ + 2 Si -^(0 jg 

positive, we have successfully transformed the partition sum to a dual representation where only real and 
positive terms appear. The proof that the sign is indeed always positive for admissible configurations of 
loops, dimers and plaquette occupation numbers will be given in the next section. 

Before we come to analyzing the sign we comment on the role of configurations with nested loops. 
We begin with the simplest case of only two loops, an outer one and an inner loop inside the outer loop. 
If the outer and the inner loop have the same orientation, the plaquette occupation numbers inside the 
inner loop will have \p{n)\ = 2. If the inner loop has the opposite orientation of the outer loop, then 
the plaquette occupation numbers inside the inner loop will be p{n) = 0, i.e., we have a hole. 


The important observation is that concerning the analysis of the sign (—+ + 2 Z^i we can 

decompose the configuration with the two nested loops into a product of two configurations without 
nested loops: One factor is the configuration of only the outer loop and the inner loop is replaced by a 
chain of dimers (see Fig. [^. This replacement of the inner loop by dimers is always possible, since a loop 
has an even number of links and every second link can be replaced by a dimer, which also constitutes an 
admissible configuration. Since the dimers all have weight factor -Fl this obviously does not introduce 
an extra sign. The second factor is the configuration with only the inner loop. If one multiplies the two 
factors, then also the plaquette occupation numbers take the correct values (see Fig. [^. 

It is obvious that the factorization trivially generalizes to configurations where we have more than 
two nested loops. The factorization allows one to analyse the sign (— 1 )^^+ + 2 Si ^(0 for the simpler 
configurations without nested loops, and in the next section we show that this sign is always -Fl. Due 
to the factorization this then implies that 




( 22 ) 
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Step 1 


Step 2 




Step 3 




Figure 3: The three steps needed for building up loops: 1) Placing the initial plaquette. 2) Attaching 
a plaquette at a single link. 3) Attaching a plaquette at three links and placing a dimer. For the 
discussion of (23) the orientation of the loop is irrelevant and we omit the arrows indicating the 
orientation. The loop can be oriented in both, the mathematically positive or the negative sense 
and the shaded areas indicate the corresponding plaquette occupation numbers +1 or —1. 


for all admissible configurations of loops, dimers and plaquette configurations. 


5 Analysis of the sign factor 


In this section we show that the sign-formula (22) holds for all configurations without nested loops, and 
thus, due to the factorization property discussed in the end of the previous section, also for all admissible 
configurations of loops, dimers and plaquettes, including those with nested loops. 

We begin the proof with showing that for a single loop I of arbitrary shape we have the following 
relation between its elements: 


2Np{l) + 2 = 4Nd{1) + 1(1) , 


(23) 


where Np(l) is the number of plaquettes inside the loop, Nd(1) the number of dimers inside the loop 
and L(l) the length of the loop, i.e., the number of links making up the loop. 


The relation (23) can be proven by building up an arbitrary loop using a finite number of iterations 


of three possible steps. For each step one can analyze the changes ANp(l), ANj:){l) and AL(l) of the 


quantities Np{l), Nd(1) and L(l) which characterize a loop. We will see that for each step (23) remains 


intact. The three steps are illustrated in Fig. For the analysis of the formula (23) the orientation of 
the loop does not play a role and we omit the arrows for the orientation of the loops in Fig. The 
shaded areas indicate plaquettes inside the loop which all have the same plaquette occupation number 
of either —1 for loops oriented in the mathematically positive sense and -|-1 for loops with negative 



















































Figure 4: Example of how a loop is built up using the three steps from Fig. 3. Above the arrows we 
show which steps were used to get to the next loop. The loop can be oriented in both, the math¬ 
ematically positive or the negative sense and correspondingly the shaded areas indicate plaquette 
occupation numbers -|-1 or —1. For simplicity we here do not show the dimers outside the loop. 


orientation. Furthermore, as an illustration in Fig. |^we give an explicit example for building up a loop 
using the three steps. 

The three steps are given by: 

1. Each loop we build starts with placing the first plaquette. This is the simplest loop one can have 
and its characteristic numbers 

Np{l) = +1 , Nd{1) = 0 , L{1) = +4 , (24) 


obviously obey (23). 

2. The loop can grow by attaching new plaquettes such that the new plaquette touches only at a 
single link of the loop. This step introduces the changes 


ANp{l) = +1 , ANd{1) = 0 , AL(0 = +2 , 


(25) 


which leave (23) intact. 

3. Finally we can create dimers in the interior of a loop by attaching a new plaquette that touches 
three links. This step comes with the changes 


ANp{l) = +1 , ANd{1) = 1 , AL{1) = -2 , 


(26) 


which also leave (23) intact. 


After the initial step 1 and a finite number of steps 2 and 3 we arrive at the loop we want to build up. 
After each step formula (23) remains intact and thus also holds for the final loop. This completes the 
proof of (23). 


For establishing the sign formula (22) we must show that the exponent on the Ihs. is even, i.e., 

1 


Np + Np + - = even . 


(27) 
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This can be done by summing the formula (23) which describes the situation for a single loop I over all 
loops in a configuration. We find 


Y^Npii) + El = 


(28) 


where we have divided (23) by a factor of 2 on both sides before summing over 1. The first sum gives 


the total plaquette occupation number Np (involve the factorization of configurations). The second sum 
is simply the number of loops Np. Rearranging the terms a little one finds 


1 


Np + Nl + -Y,L{1) = 2Y,Nd{1) + Y.L{1) = 


even 


(29) 


i i i 

where in the last step we have used that -^(0 even since all closed loops can only contain an even 


number of links. Thus (27) holds and the proof of (22) is complete. 


6 Discussion of the dual representation 

In the previous section we have shown that the sign factors are trivial. Thus we have established the 
final form of the partition sum, 




V 


E n -^|p(n)| (2\/w) 

{l,d,p} ^ 


p(n) 


(30) 


The sum J2{idp} over all admissible configurations of the dual variables, loops, dimers and 

plaquette occupation numbers. For each configuration {l,d,p) the weight factors are real and positive, 
and we have completely solved the complex action problem introduced by the topological term (as long 
as we choose rj and fj positive, i.e., as long as /3 > 1/2). 

It is interesting to discuss how the vacuum angle 6 enters the weight factors in the dual representation. 
For 0 > 0 we have rj <f] and thus yjrj/rj < 1. From (30) it is obvious that this implies that for 0 > 0, 
negative plaquette occupation numbers p(n) have a larger weight I\p{n)\{‘^Vm) (vA?7^)^^"^ and thus are 
favored. In turn this implies a larger weight for fermion loops with mathematically positive orientation. 
However, there are also pure gauge configurations that are favored: In the dual representation these are 
sheets with constant p(n) < 0 for all plaquettes, which in the conventional representation correspond to 
configurations with constant field strength, i.e., constant electric field. 

Simple bulk observables can now be obtained by taking derivatives of InZ with respect to the 
parameters fi and 6. Evaluating these derivatives also for the dual representation, the dual forms of 
the topological charge, the topological susceptibility, the plaquette expectation value and the plaquette 
susceptibility can be obtained. 

It is also possible to find the dual representations of various gauge invariant n-point functions. The 
simplest case are n-point functions of plaquettes, where one simply can introduce couplings rj{n) and 
??(n) that depend on the site index n and act as sources. After taking the corresponding derivatives of 
Z or InZ, one sets all ?7(n) and r/(n) to the desired values rj and r] and so obtains the plaquette n-point 
functions in the dual representation. 

In a similar way one can construct n-point functions of fermionic currents by multiplying all link 
variables in the fermion action ([^ by sources V/(n) and V^{n) attached to the links for forward and 
backward hops. When integrating out the fermions, the emerging loops and dimers are dressed with 
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these sources, while all other steps of the mapping to the dual variables remain unchanged. For obtaining 
the fermionic n-point functions one then evaluates the corresponding derivatives of the dual partition 
function with respect to the sources and then sets all f4(n) = Vy{n) = 1. 

We stress at this point, that the specific real and positive dualization of massless lattice QED 
presented here is possible only in two dimensions. From a technical point of view this can be seen 
from the fact that in more than two dimensions one can find simple configurations which give negative 
contributions, e.g., a configuration with a single 2x2 loop and the rest of the lattice filled with dimers. 
From a more physical point of view, the successful real and positive dualization of the lattice model 
reflects the fact that 2-dimensional fermions can be bosonized also in the continuum, i.e., they can be 
completely represented in terms of bosonic variables. 


7 Two flavor model with chemical potential 


In the conventional representation the partition sum for the two-flavor model with chemical potential is 
given by 

where il: and x si's positively and negatively charged fields respectively. The measure T’fx^x] foi' the 
second flavor x is a product over Grassmann integrals like !![?/;, V’]- The fermionic actions for the two 
flavors now contain chemical potentials and 




^ ^ \lv{n) Uy{n) ?/)(n-fz>) 

n,u 


7;y(n)e ^> 


(32) 


SxP,X,x] 


\ Y. [>(^) ^ X{n+i>) 

71 , 1 / 


Xu{n)e Uu{n)x{n+v)x{n) ■ (33) 


We stress that for the second flavor x we use the complex conjugate link variables U^{n)* = Uu{n)~^, 
since x has opposite charge. This is necessary since a theory with U(l) charges at finite density requires 
overall electric neutrality due to Gauss' law. 

We remark that for the two-flavor case considered here physics will only depend on the sum of 
the two chemical potentials +/^x Gauss’ law. However, for more than two flavors (keeping 

overall electric neutrality) physics will depend also on non-trival combinations of the individual chemical 
potentials (see, e.g., [lOj). The generalization of our dualization to more than two flavors is trivial, and 
with this possible generalization in mind we find it instructive to explicitly show the dependence on the 
individual chemical potentials in the subsequent derivation of the dualization. 

The dual representation of the two flavor model is obtained in the same way as in the one-flavor case. 
The first step is to again apply the transformation (|^ and subsequently to integrate out the fermions. 
This leads to the loop representations of the fermionic partition sums Z^[U] and Zy^lU] for the two 
flavors. The chemical potential terms give rise to extra factors of (e^^) for all forward temporal 
hops of a fermion and factors of (e“^^) for backward temporal hops. These factors cancel for 
dimers and for loops that close trivially and thus only loops I with a non-trivial winding number W{1) 
around compact time depend on the chemical potential via a factor of (i)^^ where 

Nt is the temporal lattice extent. Thus for the field the dual representation Z^[U] of the fermionic 
partition sum has the form of ( |l7| ) with an additional factor otherwise remains unchanged. 

For the field x an additional factor appears in Z^[U], and furthermore all link variables Uyiji) 

are replaced by Uy{n)~^ due to the opposite charge of x- 
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In a second step the product Z^\U] Z^[?7]e“of the two fermion partition sums and 
the Boltzmann factor for the gauge action is integrated over the gauge fields. Again the Boltzmann 
factors for the individual plaquettes are expanded using (19). Only those terms survive where all link 


variables from the fermion loops are compensated with links from the plaquettes, where for the flavor 
X all link variables along the loop are complex conjugate. The important observation is that one can 
saturate the links for the loops from the two flavors separately with the plaquettes as needed, and thus 
the result for the positivity of the overall sign proven in Section can be taken over to show that all 
configurations are real and positive also for the two-flavor case. To determine the correct weight for dual 
gauge field variables one simply adds the plaquette occupation numbers needed for saturating the link 
variables for both flavors, and this sum gives the dual plaquette variable p{n) used for the order of the 
modified Bessel function I\p(n)\ in the weight factor. 


|p(n)| 

The dual form of the partition function of the two flavor case is then given by 


'1 \ 2y 

^ = '2, 


jWJVTH'Wjl.xIVTH'd) 




p(n) 


(34) 


The sum runs over all admissible configurations of the loops and dimers for the 'j/^-field represented as 
before by I and d, over all admissible configurations of the loops and dimers for the x^held represented 
by I and d, and the corresponding admissible configurations of the plaquette occupation numbers p. For 
both flavors, i.e., for both, the l,d and the J,d variables, each site has to be either the endpoint of a 
dimer or run through by a loop. For each link the combined flux of the loops of both flavors has to 
be compensated by activated plaquettes, where the flux from the loops I that represent the oppositely 
charged field x is counted with a negative sign. 

It is obvious, that also for arbitrary values of the chemical potentials and /i^ the weights in the 
partition sum (34) are real and positive, and that in the formulation with the dual variables the sign 


problem is gone. As discussed for the case of the one-flavor model in Section]^ also in the two-flavor 
case the simplest observables are obtained as derivatives of InZ with respect to the parameters. In 
particular, here we can also consider the derivatives with respect to p^Nx and Px^^T, which give rise 
to the corresponding particle numbers for the two flavors. The form of the dependence on the chemical 


potential in (34) makes clear the interpretation of these particle numbers in terms of the dual variables: 


They are simply given by the temporal net winding numbers W{1) and W{1). As for the one-flavor 
case, n-point functions can be obtained by coupling local sources, which after taking the corresponding 
derivatives are set to 1. 


8 Possible update strategies 

We conclude our presentation of the dual representation of the massless lattice Schwinger model with 
discussing possible strategies of a Monte Carlo update. Since in the dual formulation configurations 
of the variables l,d and p (and for the two-flavor case also I and d) have to obey the admissibility 
conditions, it is not a-priori clear how to propose changes of a configuration, such that all constraints 
remain intact. In addition, the proposed changes have to be complete, such that the update is ergodic, 
i.e., all configurations can be reached with non-zero probability in a finite number of steps. We do 
not intend to discuss all technical details of possible algorithms or their performance (both these issues 
will be addressed in a forthcoming publication), but we sketch the strategies for a dual Monte Carlo 
simulation such that it is clear that the dual representation is indeed suitable for simulations without a 
sign problem. 
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We begin with noting that the dinners and the elements of the loops come with a Boltzmann factor 
of 1. For the dimers this implies, that we can shift them around for free, and only have to make sure that 
we modify them such that all fermionic constraints remain intact, i.e., all sites are either run through by 
a loop or are endpoints of dimers (for the two-flavor case this is required for loops and dimers of both 
flavors). For the loops the update is a little more involved, since changing a loop implies also a change of 
some plaquette occupation numbers p{n) because the flux along the contour of a new piece of loop has 
to be compensated by suitably occupied plaquettes. Thus for changing a loop a Metropolis acceptance 
step will be necessary, while a change that affects only dimers will always be accepted. Consequently 
the Monte Carlo strategy which we adopt here will consist of alternating updates that affect only dimers 
with local updates that change loops (and might also remove or insert a dimer). In the case of the 
two-flavor model, one has to independently do the dimer-only and the loop updates for both flavors. 
The starting configuration can be chosen as a configuration where the whole lattice is filled with dimers. 

We begin with discussing updates where the loops are modified. These updates are essentially 
patterned according to the three steps that were used in Section [5] to show the positivity of the sign 
factors and which are illustrated in Fig. The difference for using them in an update is that they can 
be applied in both directions, i.e., a loop can not only become larger (as in the positivity proof), but 
can also shrink or even vanish completely. The steps used in an update of the loops are: 

• A plaquette where two sides are occupied by two dimers is replaced by a loop and the corresponding 
plaquette variable p(n) is changed by ±1 (note that the plaquette under consideration can be nested 
inside an outer loop such that already before the change we have p{n) ^ 0. This corresponds to 
Step 1 of Fig. and of course can also be run through in both directions, i.e., we can remove an 
elementary loop around a single plaquette and replace it by two dimers. 

• A plaquette is attached to a loop along a single link and the dimer opposite that link is removed. 
The plaquette variable p{n) is changed by ±1 depending on the orientation of the loop. This 
change corresponds to Step 2 in Fig. and also can be reverted, i.e., a single-plaquette-detour 
can be removed from a loop and a dimer is placed. 

• A plaquette is attached to a loop such that it fills a one-plaquette gap in the loop. Three line 
elements of the loop are replaced by one line element and a dimer inside of the loop. The plaquette 
variable p(n) is changed by ±1 depending on the orientation of the loop. This change corresponds 
to Step 3 in Fig. and also can be reverted, i.e., a plaquette and a dimer inside a loop can be 
removed together. 

• Finally we need a step that is not depicted in Fig.[^ We can join two loops of the same orientation 
that run through opposite links of a plaquette, by deleting the two loop elements on the two links 
and inserting them on the two other links of the plaquette. Also here the plaquette variable p{n) 
is changed by ±1 depending on the orientation of the two loops. The corresponding inverse step 
cuts a loop that has a narrow section with only one plaquette into two loops. 

• For the system with chemical potential we need another step to populate the different particle 
number sectors and to explore the dependence on the chemical potential. This can, e.g., be done 
|2j by offering double lines of matter flux that wind around the compact time direction. 

All changes are accepted or rejected in a Metropolis step where the acceptance probability is computed 
from the changing weight factors I\p{n)\ of the dual gauge variables p{n) that are altered by ±1. It is 
easy to see that the steps give rise to an ergodic algorithm for the loops. 
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As outlined above, the second element of our Monte Carlo strategy is an update of only dimers while 
keeping the loops fixed in that step. The dimers have to be placed such that all sites that are not run 
through by a loop are endpoints of a dimer. As long as this constraint is obeyed the dimers can be 
placed randomly since they come with a Boltzmann weight of 1. We thus refer to our update of only 
the dimers as "random dimer insertion”. For obeying all constraints the strategy we follow is to first 
identify the dimers where we have no freedom for changing them. An example is given by the dimer 
inside the last loop of Fig. These uniquely determined dimers are frozen and will not be changed. 
For all sites where the dimers are not determined uniquely we delete the dimers attached to that site. 
Then we start to place dimers randomly connecting empty sites. After each such random placement we 
again determine whether new dimers need to be inserted that are already uniquely determined by the 
newly placed random dimers. This procedure is iterated until all dimers are placed and the constraints 
are obeyed for all sites. 

We stress that only the combination of the loop update and the random dimer insertion gives rise 
to an ergodic algorithm. Finally, to speed up decorrelation, it is possible to augment the loop and 
random dimer insertion updates with pure gauge updates, where all plaquette variables are changed by 
A G { — 1, +1}, thus inserting a sheet of constant field strength which otherwise could only be generated 
by growing a large loop all the way around both, the spatial and the temporal boundary. Concluding 
the section about possible algorithms for the dual formulation we remark again, that this discussion 
is only meant to briefly sketch the dual update strategies and a more detailed numerical analysis is in 
preparation. 

9 Concluding remarks 

In this paper we have derived a real and positive dual representation of the massless lattice Schwinger 
model with staggered fermions. The dualization solves both sign problems, the one coming from chemical 
potential and the one introduced by the topological term. The dual degrees of freedom turn out to be 
oriented loops for the fermions, filled with integer valued plaquette occupation numbers for the gauge 
fields, which can be interpreted as 2-dimensional surfaces bounded by the fermion loops. All terms in 
the partition sum come with real and positive weight factors, and we have sketched how a simulation in 
the dual representation can be implemented. To our knowledge this is the first example of a complete 
real and positive dualization of a lattice gauge theory with fermions at arbitrary values of the chemical 
potential, of the vacuum angle and of the gauge coupling (up to the trivial restriction /3 > 1/2). 

The dualization we present is for the case of massless fermions. This is an interesting fact, since so 
far most attempts for at least a partial dual representation were starting from a model at infinite fermion 
mass, i.e., the quenched case, and the dual fermion loops were obtained from hopping expansion of 
the fermion determinant, i.e., an expansion in large quark masses. For our dual representation of the 
massless Schwinger model it would be interesting to study an expansion of the dual formulation in small 
mass. For non-zero mass one finds negative contributions to the dual partition sum, but for sufFiciently 
small mass this sign problem can be expected to be mild and it could be possible to analyze the approach 
to the chiral limit in a dual representation. 

Another important lesson from dualizing this model theory is the fact that for proving the positivity 
of all weight factors the inclusion of the gauge field dynamics is an essential step. Naively one could have 
imagined that for the case of finite chemical potential, where the complex action problem comes solely 
from the fermions, one could try to find a positive representation for the fermionic partition sum alone. 
However, the construction presented here shows that the gauge field contributes exactly the necessary 
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signs to make every term in the dual partition sum positive. It can be expected, that also for other 
theories of gauge fields interacting with relativistic fermions a successful dualization strategy will have 
to consider the fermions and the gauge fields on equal footing. 
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Appendix 


In this appendix we prove the representation (19) of the Boltzmann factor for a single plaquette variable. 


For simplicity we here denote the plaquette variable as U with U G U(l). The formula we want to prove 
reads 

, _ / nfi\P 

(35) 


e-^u^-vu-^ = 


pez \ V '/ / 


where I\p\{2^/^) denotes the modified Bessel functions. The first step in the proof is to expand the 


two exponential functions on the Ihs. of (35) into their power series, 

OO OO 


-T]U ^ ^k^kjrk-k 
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(36) 


The sum over k,k £ Wq can be replaced by a sum over two new variables p G Z and p G INq. They are 
related to the original variables k and k via 


k-k = p , k + k = \p\+2p , = A: = 


Using the new summation variables we find 


e-^u^-vu-^ = 


— r>n ^ n—n 


(\/w) 


=dp|+2p 


^ {\P\ +p)'-p'- 


(37) 


(38) 


The sum over p represents a modified Bessel function (see, e.g., dl). 

(z/2)2f+’" 


In{z) = E 


{n + j)\j\ ’ 


(39) 


and replacing the corresponding sum in (38) by /|p|(2-y/r)r)) completes the proof of (35) 
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